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The possibility of an nnA bound state is investigated in the framework of pionless effective field 
theory at leading order. A system of coupled integral equations are constructed in the spin-isospin 
basis, of which numerical solutions are investigated. In particular, we make use of the limit cy¬ 
cle behavior, i.e., cyclic singularities of coupled integral equations of the system, which would be 
associated with the formation of a three-body bound state, so-called the Efimov state, in the uni¬ 
tary limit. Furthermore, we find that, when the sharp momentum cutoff introduced in the integral 
equations is taken significantly larger than the hard scale of the effective theory, the coupling of a 
three-body contact interaction becomes cyclically singular indicating the onset of Efimov-like bound 
state formation. However, the paucity of empirical information to determine the parameters of the 
theory precludes a definitive conclusion on the existence of such a bound state. As a simple test of 
the feasibility of the nnA bound system in nature, we explore the cutoff dependence of the theory, 
and uncertainties of the present study are discussed as well. 

PACS numbers: 11.10.Hi, 13.75.Ev, 21.45.-v, 21.80.-|-a. 


I. INTRODUCTION 

Recently the formation of ^n, the bound state of two 
neutrons and one A hyperon (nnA system), has been 
suggested by the experimental investigation of the Hy- 
pHI Collaboration [l|B If confirmed, this observation 
would be a crucial indication of the discovery of a new 
exotic bound state, namely, a nucleus without a proton. 
Subsequent theoretical studies , however, questioned 
the claim of Ref. [I| because of the inconsistencies of the 
putative nnA bound state with other observables such 
as the NK scattering data, hypertriton binding energy, 
and the energy gaps between the ground and first ex¬ 
cited states of and ^He. Such theoretical works 
were mostly based on the estimations employing standard 
potential model approaches, which include the KN-YjN 
mixing that play a crucial role in charge symmetry break¬ 
ing effects and determination of the spin-parity quantum 
numbers of the ground state. In fact, the seminal theo¬ 
retical work of Ref. Q already reported a long time ago 
on the non-existence of nnA bound state where model pa¬ 
rameters of a variational calculation were tuned by the 
hypertriton binding energyB 

In the present exploratory study, we investigate the 
possibility of a bound nnA system employing a pionless 
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^ The possibility of such an nnA. bound state was also suggested 
by the lattice QCD simulation of Ref. Q in the limit of flavor 
SU(3) symmetry. 

^ Other earlier theoretical studies on this system can be found, 
e.g., in Refs. [EUl- 


effective field theory (EFT) at leading order (LO) follow¬ 
ing the approach of Ref. [1^, which was used to inves¬ 
tigate the pnA system and the hypertriton bound state. 
Previously, low-ener gy E FTs were constructed by two of 
the present authors (ill. in for investigating the hyper¬ 
nuclei aaH AvHe as possible AAd and AAa bound 
states, respectively. These investigations revealed sub¬ 
stantial indications, corroborating previous claims, that 
aaH is likely to form a bound state and the bound state 
of AA^e could be an Efimov state. 

A low-energy EFT is constructed by introducing a hard 
scale Nfi that separates the relevant low-energy degrees 
of freedom from the irrelevant high-energy degrees of 
freedom which are to be “integrated out.” The advan¬ 
tage of this approach is that it provides with a model- 
independent and systematic calculational technique with 
a small number of coupling constants that embody all 
the ignorance about the short-distance dynamics. For a 
review on the details of this subject we refer the reader 
to Refs. [13, Ell and references therein. 

Throughout this work, we are dealing with a three- 
body system which, if bound, is likely to have a binding 
energy much smaller than the pion massB Therefore, we 
can choose the pion mass as the hard scale, i.e., Kh ^ 
so that the pions are integrated out and are not 
explicitly introduced in the theory. We may additionally 
regard the mass difference between the A and E hyperons, 
i.e., 5m = m-^ — ~ 80 MeV, as of the same order as 

the hard scale Kh- Then our effective Lagrangian can 
be written explicitly in terms of the neutron and A fields 


® The estimated binding energy of the putative nnA bound state 
in Ref. [J is at most a few MeV. 
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along with their interactions described by the contact 
terms, which will be determined in a phenomenological 
way. 

In addition, we will make use of the cyclic singularities 
that arise in the solutions for the coupled integral equa¬ 
tions in the asymptotic limit d- Such singularities are 
renormalized by introducing a suitably large momentum 
cutoff Ac (Ac > Ah) in the loop integrations at the cost 
of introducing three-body counter terms at LO. Conse¬ 
quently, in order to absorb this cutoff dependence, the 
corresponding three-body coupling may exhibit a cyclic 
renormalization group (RG) evolution termed as the limit 
cycle d- The cyclic singularities are associated with the 
occurrence of bound states, known as the Efimov states, 
in the resonant/unitary limit d- In our analysis, we 
vary the magnitude of the cutoff within a reasonable 
range to investigate its sensitivity to the formation of 
bound states. This is a simple test one can perform for 
checking the feasibility of the nnA system as a three-body 
bound system. The main purpose of the present inves¬ 
tigation is, thus, to explore the putative bound state of 
the nnA system in the context of a modern EFT method 
mentioned above. 

This paper is organized as follows. In Sec. im the ef¬ 
fective Lagrangian at LO for the system is introduced. 
Then the renormalized dressed two-body propagators are 
defined and the three-body coupled integral equations 
are derived in Sec. uni Before numerically solving the 
coupled integral equations, we obtain an analytical ex¬ 
pression of a scale-invariant equation needed to examine 
the limit cycle behavior in the asymptotic limit by as¬ 
suming that the neutron and A hyperon have the same 
mass. Section m is devoted to numerical solutions of 
the coupled integral equations for the nnA system with 
physical baryon masses, and we test if the three-body 
contact interaction exhibits cyclic singularities, typically 
associated with the Efimov states, even at low and in¬ 
termediate momenta being away from the unitary limit. 
We further investigate the possibility of bound state for¬ 
mation in the absence of the contact interaction, when 
the sharp cutoff Ac is chosen to be significantly larger 

I 


than Ah- Section El summarizes the present work with 
possible implications of our results and the uncertainties 
involved in the present approach. 


II. EFFECTIVE LAGRANGIAN 


The relevant non-relativistic effective Lagrangian for 
the nnA system at LO consistent with parity, charge con¬ 
jugation symmetry, time-reversal invariance, and small- 
velocity Lorentz transformation reads 


£ — -f £a 






(nA) 


■'t(nA) 


c 


3-body 


, ( 1 ) 


where the elementary fields of our EFT are the neutron 
field Bn and the A hyperon held Ba of which one-body 
Lagrangian are represented by £„ and £^, respectively. 
The S'-wave two-body Lagrangian for the spin singlet nn 
channel, spin singlet nA channel, and spin triplet nA 
channel are respectively represented by £s(„„), 
and Ct(nA) ■ The composite dibaryon helds are introduced 
and denoted by S(„„) for the spin singlet nn, and S(„a) 
and t(„A) for the spin singlet and spin triplet nA systems, 
respectively. The three-body interaction Lagrangian is 
represented by /la-body 

The one-body Lagrangian £„ and of Eq. ([T]) are 
given as [H, [13 
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( 3 ) 


where is a velocity four-vector chosen as = (1,0). 
The neutron and A hyperon masses are given by and 
TO^, respectively. The ellipses denote higher order terms, 
which are not required at the accuracy of the present 
analysis. 

The S-wave two-body interactions are written in terms 
of the dibaryon fields as 


p _ f 

>~-s{nn) — ^s(nn) ^(nn) 


_ t 
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H.c. 


( 4 ) 
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( 6 ) 


where <^s(nA)^ and are sign factors, and 

Ag(„A)) and A((„a) are the respective mass dif¬ 
ferences between the corresponding dibaryon and its con¬ 
stituent elementary particles. The coupling constants are 


denoted by and which will be de¬ 

termined by the S'-wave effective range parameters such 
as the scattering lengths and effective ranges [20l - [2^ as 
will be discussed in the next section. The spin projection 
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operators introduced in the above Lagrangian are defined 
as 


P, 


:'so) 


(nn) 


--a 

2^2 J ^(nh) 


= - 
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where the difference in the factors in the definitions of 
the spin singlet projection operators, and 

arises from the existence of two identical particles in the 
nn channel. 

As will be demonstrated in the next section, the so¬ 
lutions for the coupled integral equations in the spin 
doublet channel of the nnA system exhibit cyclic sin¬ 
gularities in the asymptotic limit. To renormalize the 
cyclic singularities, it mandates the inclusion of three- 
body counter terms already at LO, where the singular¬ 
ities associated with the short range part of the one- 
neutron/A-exchange interactions are needed to be can¬ 
celed out. In three-nucleon systems, the general expres¬ 
sion for three-nucleon counter term is constructed in the 
Wigner SU(4)-symmetric limit [l^, [2^. Moreover, in the 
case of the hypertriton, i.e., the spin doublet npA chan¬ 
nel, the general expression for the three-body counter 
term is obtained in the limit where the mass difference 
between the nucleon and the A hyperon is ignored [la. 
Following this approach, as will be shown in the next sec¬ 
tion, one may construct an analogous expression for the 
counter term using the projection operator for the spe¬ 
cific diagonal mode required to renormalize the asymp¬ 
totic cyclic singular behavior in the same mass limit. 
However, with the physical value of the mass difference 
taken into account, a more systematic way of 
analytical determination of the counter term is challeng¬ 
ing. Furthermore, we will make use of the RG evolution 
of the bound states that exhibits the periodic occurrence 
of critical cutoff values, corresponding to the vanishing 
points of the three-body contact interaction, in studying 
the limit cycle behavior of this system at threshold, i.e., 
with zero three-body binding energy. 

In our present investigation, therefore, we introduce 
the three-body counter-term Lagrangian for the spin- 
doublet channel as 


Tla-body — g™A2/t(„A) ^2 *(nA)i 
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where g{Ac) is a cutoff-dependent coupling constant. In 
general, the 3-body contact interaction of Eq. con¬ 
tains other terms that involve the nn dibaryon and A hy¬ 
peron fields as well as the spin singlet nA dibaryon and 
neutron fields with different coupling strengths. How¬ 
ever, as will be discussed later, the paucity of data on 
the nnA system at the present stage does not allow the 
estimation of these couplings. Thus, the term shown in 
Eq. ([5]) will adequately serve our purpose in this explo¬ 
rative study where we hope to capture certain universal 


features of such bound systems even without resorting 
to a more sophisticated EFT analysis. In principle, the 
above coefficient g{Ac) should be renormalized as a func¬ 
tion of the same cutoff parameter Ac that is introduced in 
the coupled integral equations for the three-body system 
and thereby fixed by the experimental/empirical data on 
a three-body observable of the nnA system such as the 
three-body binding energy. 


III. AMPLITUDES AND COUPLED INTEGRAL 
EQUATIONS 

Throughout this study we employ the power counting 
rules suggested by Kaplan, Savage, and Wise (KSW) for 
the two-body sector [^[^. At LO we consider the S- 
wave two-body nn and nA interactions that leave us with 
the three channels in the two-body sector, namely, the 
nn in ^Sq state, nA in ^Sq state, and nA in state. 
The KSW rules require the “bubble” diagrams in the 
two-body propagators of nn and nA to be resummed to 
infinite order, while the propagators are renormalized at 
LO using a single parameter, i.e., the respective S'-wave 
scattering length. 

For the three-body sector, we generally follow the 
prescription suggested by Bedaque, Hammer, and 
van Kolck Namely, a three-body contact interac¬ 
tion, which would otherwise be naively considered as a 
subleading contribution, should be promoted to LO. This 
is because of the non-analytic ultra-violet (UV) enhance¬ 
ments from the cutoff dependence which appears when¬ 
ever the three-body system exhibits the limit cycle be¬ 
havior. Since there is no empirical information to fix 
the strength of the three-body interaction, however, the 
three-body contact interaction is taken into account only 
when the limit-cycle behavior is explored, while we as¬ 
sume a vanishing binding energy of the nnA system 0 
When we explore the nnA binding energy, on the other 
hand, we exclude the contributions of the three-body con¬ 
tact term by focusing on the role of the two-body inter¬ 
actions. In other words, instead of making a definite 
prediction on the binding energy of the nnA system, we 
try to reveal some general features of the system. 

For the two-body sector, the Feynman diagrams for 
the dressed dibaryon field in the nn(^S'o) channel are 
shown in Fig.[Tl which leads to the renormalized dressed 
dibaryon propagator at LO as 

Ds(nu)[%,(l) = -2 / , „ ’ 

'/*(„„) n -- \J zQ - - ze 

(9) 

where is the scattering length of the neutron- 

neutron scattering in the channel, and Qq and q 


^ Therefore, this corresponds to the case when the energy of the 
system equals to its threshold energy. 
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FIG. 1. Diagrams for the dressed dibaryon propagator for the spin singlet nn channel. The single and double lines denote the 
neutron and dibaryon fields, respectively. 


are generic off-shell energy and three-momentum, respec¬ 
tively. The loop diagrams are calculated using dimen¬ 
sional regularization with the power divergence subtrac¬ 
tion scheme [2^ [2^ which introduces the subtraction 
scale parameter fi. The coupling constants in the La- 
grangian are renormalized by using the S'-wave effective 
range parameters, namely, the scattering length and 


1 


the effective range r„„. Furthermore, the analyses in 
Refs. (20l - [^ yield the sign factor = — 1 and 
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47rA, 


(nn) 


^nu rn^Vsinn) 


Stt 
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Similarly, for the nA channels. Fig. [5] displays the 
Feynman diagrams for the dressed nA propagators whose 
renormalized expressions are given byl£| 


^s{nA) {Qo 1 

for the spin singlet channel and 


27r 


ys{nA)^^{nA) 


1 

^s(nA) 


Dt(nA){%,q) 


27r 


^t(nA)l^(nA) 


1 

^i(nA) 


for the spin triplet channel, respectively, where cLs(^nA) 
(«t(nA)) is the S'-wave scattering length of nA scatter¬ 
ing in the spin singlet (triplet) channel, and ii(„A) is the 
reduced mass of the nA system. Again the coupling con¬ 
stants of the two-body Lagrangian are renormalized by 
the S-wave effective range parameters and, as detailed in 


Refs. 


, we have a 
1 

^s,t(nA) 
^s,i(nA) 


s,t(nA) 


= — 1 and 


27rA, 


,t(nA) 


^(nA)^s,t(nA) 

27r 

2 2 * 
^(nA) ^s,t(nA) 


(13) 


For the three-body part, we derive a set of coupled in¬ 
tegral equations for the S-wave scattering of the neutron 
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■ l/-2M(nA) 1 Qo - 
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and the composite dibaryon triplet state t(nA) in the mo¬ 
mentum space. There are two possible allowed total-spin 
channels, namely, spin-3/2 and spin-1/2. Since the inte¬ 
gral equation for the spin-3/2 channel does not exhibit a 
limit cycle0 we shall henceforth consider only the spin- 
1/2 channel as demonstrated by the diagrams in Fig. |31 
The coupled integral equations are expressed in terms 
of the three half-off-shell amplitudes, a{p',p), b{p',p), 
and c{p',p), where a{p',p) is the amplitude of elastic 
nt{nA) scattering, b{p',p) is that of inelastic nt{nA) to 
ns(nA) scattering, and c{p'p) is that of inelastic nt{nA) 
to As(nn) scattering, with p' (p) being the relative off- 
shell (on-shell) momentum of the final (initial) two-body 
system. Explicitly, we have 


dlfK^,^{E-p',l) - 

1 

SCnA) 


b{l,p) 



2M(nA)-® 


a{l,p) 





E 


^ The ie prescription in the dressed propagators is understood in 
the same way as in Eq. (0 for Ds(„„)(gp,i 5 f). 


In the spin-3/2 channel, one obtains a single integral equation of 
which numerical solution does not yield the limit cycle behavior. 
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+ Kl>=C)= + ■ 

FIG. 2. Diagrams for the dressed dibaryon propagator for the spin singlet and triplet nA channels. The double line denotes 
the dibaryon field and the thin (thick) line the neutron (A hyperon) field. 





FIG. 3. Diagrammatic representation of the coupled integral equations for the neutron and spin triplet nA dibaryon elastic 
scattering in the spin-doublet channel without three-body contact interaction. Blobs with the horizontal lines denote the elastic 
amplitudes for the neutron and spin triplet nA dibaryon [nt(nA)] channel, while those with the vertical and crossed lines denote 
the two inelastic channels, namely, nt(nA) to neutron and spin singlet nA dibaryon, and A hyperon and the spin singlet nn 
dibaryon channels. See the captions of Figs. [T] and [2] as well. 




^/3 


®nn 


c{l,p) 


2fi 


-P — m„E 


(nn)A 


m. 


b{p',p) = -7r^Kyt{nK)Kl^a){E]p',p) - . 

^ P‘(nA) JO 


dlPK^a)iE-,p',l)- 


a{l,p) 


“■t(nA) 


-^P-2p, ..E 


a(nA) 


+ 


I m 


A 


2 ^ PinA) Jo 


dlPK^a){E-p',l)- 


b{l,p) 


^s(nA) 


Ei^P-2p, ,,E 

^n(nA) ^(nA) 


+ ^ dlPK^t2){E;p',l) — 

^ Jo — 


c{l,p) 


2^ 


-P — m„E 


(nn)A 


c{p,p) = \l -jm^yt(nA)K(bi){E-,p',p) - 


3 1 


2 P-{nA) Jo 


dlPK(bi){E-p',l) 


i{l,p) 


(14) 


(15) 


“t(nA) Y '""("A) 


+ 


1 




nA) *^0 


dlPK(bi){E-p',l)- 


-‘s(nA) y ^n(nA) 


K^.p) _ 

tE^P-2p, .,E 

. . .. ^(nA) 


(16) 


where A(a)(i?;p',p) is the one-A-exchange interaction kernel, and Kp,i){E\p' ,p) and K(i, 2 )iE;p',p) are the two possible 
one-neutron-exchange interaction kernels with E being the total center-of-momentum energy of the three-body system. 
These interaction kernels are written as 
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K(bi){E;p',p) 


K[b2){E-,p\p) 



where Pn{nA) and pi{nn)A are reduced masses defined as 
A<n(nA) = + Tnj^)l{2m„ + m^) and p(nn)A = 

2m„TO^/(2m„ + m^). As mentioned earlier, the sharp 
momentum cutoff Ac is introduced in the above integral 
equations assuming E,p ^ l/Orm) ^/^s,t(nA) ^ P' ^ ^c- 
While the inhomogeneous parts of the integral equa¬ 
tions only set the overall low-energy scale of the prob¬ 
lem in the asymptotic limit of the off-shell momenta 
p' and I, i.e., p',l ^ Ac, the behavior of the solutions 
are completely determined by the homogeneous parts. 
This, however, leads to ambiguities in the solution when 
the equations are numerically solved in the UV limit 
(Ac —>■ oo) without the three-body contact interaction. 
As evident in the analysis below, we need to introduce 
a three-body counter term only for one of the three 
elastic scattering modes in the diagonal basis associated 
with the asymptotic limit cycle behavior. Nevertheless, 
when the integral equations of Eqs. (O-dini) are numeri¬ 
cally solved in the spin-isospin basis for intermediate mo¬ 


p'^ +p^ + p'p -m^E 


.p/2 _|_ p 2 _ p/p _ rn^E j ’ 

p^ + p'p — m„E\ 
-p^ — p'p — m^E j ’ 


2/j. 


(nA) 


2/J. 
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(18) 

(19) 


menta, one requires three-body contact interaction terms 
with several different unknown couplings. In this prelim¬ 
inary analysis, this is represented by Eq. ®, and we only 
need to modify the one-A-exchange kernel in the first two 
terms in Eq. da as 


K^,^{E-p',l)^K^,^{E-p',l) 


g(Ac) 

A? 


( 20 ) 


Before discussing the numerical solutions of the cou¬ 
pled integral equations of Eqs. (HI- (ll^. let us consider 
the approximation suggested in Ref. [lOj for investigat¬ 
ing the asymptotic nature of the solutions excluding the 
three-body interaction. Namely, since the mass differ¬ 
ence between the A hyperon and neutron is small, i.e., 
6 = (m^ — TO„)/(m^ -I- TO„) ~ 0.1, the corrections due to 
small 5 to the integral equations may be ignored in the 
asymptotic limit. Then we have 
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where 


^ 2p'l^''[p'^+E-p'lJ ' 


( 22 ) 


In the asymptotic limit of the respective amplitudes, the 
dependence on the on-shell momentum p is implicitly un¬ 
derstood. The integral equations can then be diagonal¬ 
ized to obtain generic homogeneous eigenvalue equations 
as 

A / /J dl fp''^ + l'^ + p'l\ . 

where A„ with n = 1,2,3 are the eigenvalues of the above 
3x3 matrix obtained as 


(24) 


■'"= 75 ”n^* 


Numerically, we have Ai ~ 3.54, A 2 ~ 0.80, and A 3 


—4.34. The amplitudes A„(p',p) are given by 
An(p') = - v^An - A„^ a(p') 
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with an arbitrary normalization. 

The expression in Eq. (1^51) has no scale dependence. 
The scale invariance in the asymptotic limit, p', I Ac ^ 
p, suggests that the S'-wave projected amplitudes A„(p') 
must exhibit a power-law behavior as [27| 


Mp')^p'^-\ 

Then, by the Mellin transformation, Eq. 
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An sin (f s) 
s cos (^s) ’ 


(27) 


where the solutions, in general, are complex-number 
functions and always come in pairs due to the additional 
inversion symmetry of A{p') —>■ A{\/p') in Eq. (0^ . 
which is equivalent to a; —>■ 1/x in Eq. (I27|) . Thus, 
there can be up to a quadruplet set of solutions, i.e., 
{±s, ±s*}, to the algebraic part of Eq. ([27|) . all of which 
may not solve the integral equation. Of these, the 
physically acceptable solutions correspond only to those 
amplitudes which converge as p' —>■ oo and for which 
the Mellin transformation exists This scale- 

invariant equation then gives the condition to examine 
the limit cycle behavior of the spin-doublet nnA sys¬ 
tem in the asymptotic limit. It is found that whenever 
the value of A„ exceeds the critical value Ac, i.e., when 
A„>Ac = 6/7r«1.91 [13, the solutions become imagi¬ 
nary and a limit cycle appears in the system. Thus, only 
the first mode with Ai ~ 3.54 should correspond to a 
limit cycle behavior with an imaginary solution, s = isQ 
with Sg = 0.80339 • • • . In addition, for the purpose of 
renormalization it may be convenient to project out this 
mode from the other diagonal modes that do not exhibit 
the limit cycle. To this end, the expression of the counter 
term in the limit, where the mass difference between neu¬ 
tron and A hyperon is ignored, can be constructed by 
considering the projection operator for the first mode, 
i.e., {M — \ 2 l){M — X 3 I), where M is the 3x3 matrix 
appearing in Eq. m and I is the 3x3 unit matrix. 

In the following section we investigate the numerical 
solution of the coupled integral equations of Eqs. 0 - 
(HH) valid for the non-asymptotic momentum range p < 
p' <C Ac with the physical A-n mass difference. Our 
analysis will reveal that the characteristic cyclic behavior 
reminiscent of the asymptotic limit cycle clearly survives 
at low and intermediate momentum range. The solutions 
exhibit the strong cutoff dependence and are very likely 
sensitive to the short-distance dynamics. This is similar 
to what was seen, e.g^in the case of the triton as well 
as the hypertriton channels. In contrast, the other 
two modes have well-behaved UV stable solutions that 
do not need a leading-order three-body renormalization. 
These modes do not exhibit the limit cycle and insensitive 
to the short-distance physics as seen, e.g., in the case of 
the spin quartet channel of nd scattering [ 21 , 133 . 


IV. NUMERICAL RESULTS 

In this section, we present our numerical results and 
investigate the singularities of the coupled integral equa¬ 
tions for the nnA system. Our EET has four param¬ 
eters at LO to be determined in a phenomenological 
way. They are the nn and nA S'-wave scattering lengths, 
namely, a„„, and and the three-body cou¬ 

pling 5 (Ac) with the cutoff Ac. In the present study, 
we determine these parameters as follows. Eirst, we use 



FIG. 4. (Color online) Strength of the coupling g{Ac) of the 
three-body contact interaction as a function of the cutoff Ac, 
where zero binding energy {B = 0) of the three-body nnA 
system is maintained. 


the standard experimentally extracted value of the nn 
scattering length, a„„ ~ —18.5 fm [sij- Since there is 
no empirical data for nA scattering, however, we quote 
the theoretically calculated values of and re¬ 
ported in Ref. [13 , which are obtained based on the 
potential constructed in the chiral effective theory up 
to next-to-leading order (NLO). Then, invoking isospin 
symmetry leads to the estimates: as(„A) — “2.90 fm and 
~ —1.60 fm. Einally, there is no experimental 
or empirical information to fix the strength of the con¬ 
tact interaction g{Ac). Thus, we investigate the limit 
cycle behavior tuning the energy of the nnA system to 
the threshold value, i.e., to vanishing three-body binding 
energy B = 0. 

In Eig. m we plot the strength of the three-body con¬ 
tact interaction g{Ac) as a function of the cutoff Ac, 
which gives zero binding energy of the three-body system. 
The periodic RG evolution of g{Ac) clearly indicates that 
the short-range part of the one-baryon-exchange interac¬ 
tions become singular, and thus the contact interaction is 
needed at LO for the renormalization of the singularities. 
In addition, one finds that the interactions without the 
three-body force become stronger for larger values of the 
cutoff Ac. However, in the region of small Ac Ah, the 
two-body interactions alone cannot generate attractive 
forces strong enough to form a bound state. This is man¬ 
ifest through the large negative values of (/(Ac) required 
to supplement the deficit in attractive forces for bound 
state formation. The lowest critical value of Ac, where 
(/(Ac) changes its sign in the above scenario, is found to 
be at around 1.5 GeV. Consequently, if Ac > 1.5 GeV, 
an nnA bound state can be generated from two-body dy¬ 
namics alone without the necessity of introducing three- 
body terms. On further increasing Ac, the second critical 
cutoff appears at Ac ~ 80 GeV, which indicates the onset 







FIG. 5. (Color online) Binding energy of the nnA three-body 
system as a function of the cutoff Ac without the three-body 
contact interaction. 


of a second bound state. 

Based on the universality of the Efimov-like states, 
the occurrence of periodic critical points is expected as 
g{An+i) = g{Ai) = 0 for A„+i = Ai exp (riTr/so)- The 
result shown in Fig. |4] reveals that Ai ~ 1.54 GeV and 
A 2 ~ 77.8 GeV, which leads to Sq = 7r/ln(A2/Ai) Ri 
0.801. This is surprisingly in good agreement with the 
value of Sq obtained in the previous section, suggesting 
an inherent universal behavior of the nnA system gov¬ 
erned by the asymptotic limit cycle behavior that is not 
‘washed away’ even for non-asymptotic momenta without 
considering the degenerate mass approximation. 

We now turn to the cutoff-dependence of the binding 
energy of the nnA system without the contact term, i.e., 
by setting g{Ac) =0. In Fig. [5] we display the result 
for the the three-body binding energy B without includ¬ 
ing the contact interaction. Glearly, the first three-body 
bound state appears at Ac ^ 1.5 GeV, and the second 
one at Ac 80 GeV. This shows the periodic nature of 
the binding energies as a function of Ac. As seen in the 
figure, the binding energies of the shallow bound states, 
as they are formed, progressively increase from zero with 
increasing Ac, and eventually becoming very large, asso¬ 
ciated with deeply bound states. 

V. SUMMARY AND DISCUSSION 

In the present study we investigated the possibility of a 
bound nnA system in the pionless EFT at LO by explor¬ 
ing the structure of the coupled integral equations de¬ 
scribing the system. Because of the limited information 
to fix the parameters of EFT, we are unable to arrive at 
a definitive conclusion on the formation of a bound nnA 
state. However, by numerically solving the coupled inte¬ 
gral equations including the three-body interaction with 
the binding energy fixed to a certain three-body thresh¬ 


old value, e.g., H = 0 as used in this study, the coupling 
g{Ac) was found to develop cyclic singularities, which is 
a characteristic of Efimov-like bound states. In addition, 
we studied the role of short range two-body mechanisms 
involving the one-A and one-neutron exchange interac¬ 
tions by setting g{Ac) = 0. Our analysis reveals that 
bound states appear without the requirement of the con¬ 
tact interaction only if the cutoff was chosen unnaturally 
large, Ac > 1.5 GeV ;#> Ah, which indicates the non¬ 
trivial role of short-range two-body mechanisms that are 
beyond the scope of this work. 

As mentioned before, the predictability of our ap¬ 
proach relies on the knowledge of the values of four low- 
energy parameters. In the present study, we use the em¬ 
pirically known value of the nn scattering length, i.e., 
®nn ~ —18.5 fm [3l|. Since there is no data for the nA 
scattering lengths, we invoke isospin symmetry to deter¬ 
mine their values from the pA scattering lengths. How¬ 
ever, the poor statistics of experimental data do not al¬ 
low precise partial wave analyses and the estimated scat¬ 
tering lengths have inevitably large uncertainties. For 
example, in Ref. [s^, the ranges of the pA scattering 
lengths are given as as(nA) — (—1.85 ^ —2.78) fm and 
~ (—1.04 ~ —1.90) fm depending on the model 
potential. In the present study we used the NLO chiral 
EFT results reported in Ref. [^: as(nA) = —2.90 fm and 
o-t{nA) = —1.60 fm. However, the previous version of the 
chiral EFT calculation at LO reported significantly small 
values for the scattering lengths, — —1-91 and 

~ —1.23 fm [s^. (See also Ref. [S^.) Thus, the 
currently estimated values of the pA scattering lengths 
are beset with certain uncertainties which evidently prop¬ 
agate into our analysisQ For instance, if we use the 
values of Ref. [S^, the first bound state of the system 
appears at Ac ~ 2 GeV, which should be compared with 
Ac ~ 1.5 GeV obtained with the values of Ref. [S^. How¬ 
ever, this does not have any significant effect on the cyclic 
nature of the bound system. 

In contrast with the case of the scattering lengths, 
there is no a priori way to determine the strength of the 
three-body contact interaction without an essential in¬ 
formation on a measured three-body observable. Thus, 
owing to the absence of enough information to fix the 
parameters of the EFT developed in this work, our cur¬ 
rent results are not conclusive on the existence of the 
putative nnA bound state. Since we treat the effective 
NA interaction as point-like, short range mechanisms 
such as the two-pion exchange interactions appearing as 
box diagrams with NA-NH mixing, or other meson ex¬ 
changes such as the rj and K are to be implicitly taken 
into account by the contact interactions. Thus, without 
a three-body contact interaction, a reasonable value of 


^ From Refs. [3^. l35ll . we find that the NA scattering lengths con¬ 
verge to values around —2.9 to —2.5 fm for the singlet case and 
— 1.7 to —1.4 fm for the triplet case. 
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the cutoff would be commensurate with the energy scale 
of the two-pion/Ty-meson/Jf-meson exchanges of the NA 
interaction, which leads to Ac 300 — 500 MeV. This 
value is consistent with the values of Ac found in sim¬ 
ilar context in the literature. For example, when the 
triton and aaKb are treated as three-nucleon and AAa 
bound systems, respectively, the cutoff values are ob¬ 
tained as Ac ~ 380 MeV for the case of the triton 
and Ac ~ 400 — 570 MeV for the case of aaHc T^ . 
depending on the adopted value of the ^-wave AA scat¬ 
tering length with the three-body contact interactions 
turned off. This contrasts with our case as we indeed 
need a very large cutoff to a form bound state with¬ 
out the three-body contact term. Therefore, the short- 
range two-body attractive mechanisms discussed above 
are unlikely to generate enough attractions in the one-A- 
exchange or one-neutron-exchange interaction channels 
that could lead to nnA bound states. This ultimately 
would make it difficult to form an nnA state in nature. 
However, further rigorous investigations, both theoreti¬ 
cal and experimental, are necessary to pin down to true 
character of the nnA system. 
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